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Abstract: We solve a random two-matrix model with two real asymmetric matrices 
whose primary purpose is to describe certain aspects of quantum chromodynamics with 
two colours and dynamical fermions at nonzero quark chemical potential fi. In this sym- 
metry class the determinant of the Dirac operator is real but not necessarily positive. 
Despite this sign problem the unquenched matrix model remains completely solvable and 
provides detailed predictions for the Dirac operator spectrum in two different physical sce- 
narios/limits: (i) the e-regime of chiral perturbation theory at small fi, where //^ multiplied 
by the volume remains fixed in the infinite-volume limit and (ii) the high-density regime 
where a BCS gap is formed and /x is unsealed. We give explicit examples for the complex, 
real, and imaginary eigenvalue densities including Nf = 2 non-degenerate flavours. Whilst 
the limit of two degenerate masses has no sign problem and can be tested with standard 
lattice techniques, we analyse the severity of the sign problem for non-degenerate masses 
as a function of the mass split and of /i. 

On the mathematical side our new results include an analytical formula for the spectral 
density of real Wishart eigenvalues in the limit (i) of weak non-Hermiticity, thus completing 
the previous solution of the corresponding quenched model of two real asymmetric Wishart 
matrices. 

Keywords: Spontaneous symmetry breaking, matrix models, chiral Lagrangians, lattice 
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1 Introduction 

The phase diagram of quantum chromodynamics (QCD) as a function of temperature T 
and quark chemical potential fi is of great importance for many phenomenological questions 
in, e.g., heavy-ion collisions, the physics of neutron stars, and the early universe. Clearly, 



- 1 - 



nonperturbative methods are needed to study the phase diagram. Lattice QCD simula- 
tions have obtained solid quantitative results at nonzero T, see, e.g., [1], but what happens 
at nonzero fi is much less well established because lattice simulations are hindered by the 
fermion sign problem, see [2] for a review. Some exact results could be derived for asymp- 
totically large fi, see, e.g., [3] for a review, but otherwise studies of the ^-dependence of 
the QCD phase diagram had to rely on model approaches. Examples include the instanton 
liquid model [4], the PNJL model [5], and universality arguments [6]. However, there are 
QCD-like theories in which the fermion sign problem can be avoided for certain choices of 
the parameters, and which can therefore be simulated on the lattice. An example is QCD 
with two colours, where the single-flavour fermion determinant is real but not necessarily 
nonnegative so that the sign problem is absent for an even number of quark flavours whose 
masses are pairwise degenerate. Such theories are important testing grounds for effective 
theories or models since a direct comparison of analytical results with lattice data is pos- 
sible. If the parameters of the theory are detuned from their "sign- free" values (e.g., if the 
degeneracy of the quark masses is lifted in two-colour QCD), a sign problem develops. As 
long as the sign problem is weak enough, one can still obtain some information from the 
lattice data, e.g., by reweighting, and therefore test the effective theory or model also in 
the presence of the sign problem. In fact, lattice studies of this kind have already been 
done with two-colour adjoint staggered fermions [7, 8], and we will comment on them at 
the end of section 4.2.2. 

In this paper we will focus on two-colour QCD with Nf quark flavours in the funda- 
mental representation of the gauge group at nonzero fi.^ This theory has been investigated 
in the past in great detail for both small and large /i. For small fx, it was studied using an 
effective chiral Lagrangian, see, e.g., [9], and in lattice simulations, see, e.g., [10]. One of 
the main results from these studies is a transition from a phase with a chiral condensate to 
a phase with a diquark condensate at a critical value of fi equal to one half of the Nambu- 
Goldstone boson mass. For asymptotically large fi, the physics is completely different. The 
ground state is believed to be a BCS superfluid of diquark pairs, for which the BCS gap 
has been computed in a weak-coupling approach in [11, 12]. The effective Lagrangian for 
this case has been derived in [13]. 

A very useful analytical approach to QCD and QCD-like theories is random matrix 
theory (RMT), which can be used either as an effective theory that yields exact results 
in well-defined limits, or as a schematic model to obtain qualitative results. RMT has 
applications in many areas of physics, see [14-21] for reviews of some of these areas. The 
RMT approach to QCD was pioneered in [22], and reviews can be found in [17, 23, 24]. 
Random matrix models for QCD at nonzero have been proposed using a single random 
matrix [25] or two random matrices [26]. Both models describe the same physics, but the 
two-matrix model turns out to be easier to handle mathematically [27], due to the fact 
that a representation in terms of complex eigenvalues exists. Random two-matrix models 
were proposed for adjoint QCD at nonzero fj, [28, 29] and most recently for two-colour 

^Most of our results apply not only to two-colour QCD but also to any QCD-like theory with pseudo-real 
quarks. 
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QCD at nonzero fi [30], completing the set of all three chiral symmetry breaking patterns 
[31]. Each RMT yields exact results in the so-called microscopic regime, corresponding to 
the lowest order in the e-expansion [32] of the chiral effective theory. For a comparison to 
lattice simulations (when these are possible) we refer to [24] . 

For the case of two-colour QCD, which is the topic of this paper, the random two- 
matrix model is explicitly given in eq. (2.1) below. It implicitly contains two parameters, 
see eq. (2.4). One parameter relates the random-matrix quark masses to the physical quark 
masses and can be identified with the low-energy constant G (in the notation of [9]) in the 
chiral Lagrangian. The other parameter relates the random-matrix chemical potential to 
the physical chemical potential and can be identified with another low-energy constant, F. 
Computing analytical results from RMT and matching them to lattice data then allows one 
to determine these two low-energy constants. It was realized in [33] that the model for two 
colours proposed in [30] also describes the completely different physics at large fi, provided 
that the random-matrix parameter corresponding to fi is set to a particular value which 
corresponds to maximal non-Hermiticity. The remaining parameter can then be identified 
with the BCS gap A. By matching the RMT results to lattice data one could in principle 
determine A, but it is currently unclear what computational resources are needed to go to 
large enough values of // on the lattice. 

The random two-matrix model proposed in [30] was solved in [34] for Nf = only. 
The eigenvalues of the pseudo-real Dirac operator are either purely real, purely imaginary, 
or come in complex conjugate pairs. The spectral correlation functions of these eigenvalues 
were computed in the quenched case (i.e., for Nj = 0) for finite N, where N is related to the 
dimension of the random matrix and assumed to be proportional to the four- volume V4. For 
N ^ 00 two different mathematical limits were considered, that of weak non-Hermiticity 
[35, 36] (corresponding to /j^'^N finite) and that of strong non-Hermiticity (corresponding 
to nonzero /i in the N ^ 00 limit). In the strong non-Hermiticity limit, the microscopic 
spectral density was obtained for the complex, real, and imaginary eigenvalues. In the 
weak non-Hermiticity limit, the microscopic spectral density could only be obtained for the 
complex eigenvalues. In this paper, we generalise the results of [34] to nonzero A'^^ using the 
building block for finite N from [37], and we also obtain the microscopic spectral densities 
of the real and imaginary eigenvalues in the weak non-Hermiticity limit. Furthermore, we 
analyse the sign problem that arises when the quark masses are detuned for even A'^^ or 
when fi is increased for odd NjP It has been argued for three colours that the sign problem 
and the resulting strong oscillations of the spectral density are directly responsible for the 
breaking of chiral symmetry in the unquenched theory [39, 40]. 

In the remainder of this paper, we will sometimes refer to the two mathematical limits 
of weak and strong non-Hermiticity as "low density" and "high density" , respectively, with 
the understanding that by "low density" we mean the physical regime in which fi'^F'^V/i is 
finite and by "high density" we mean the physical regime in which a BCS gap is formed. 
Our RMT results can be directly applied to lattice simulations of two-colour QCD in both 

■^The possibility of unequal chemical potentials as yet another source of the sign problem is interesting, 
see e.g. [38], but will not be considered in this paper. 
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regimes. 

The structure of this paper is as foUows. In section 2 we define the partition function in 
RMT as well as in the two limiting cases of low and high density. New Pfaffian expressions 
for the latter are given. Section 3 first summarises the structure of finite-A^ results for 
all spectral densities, expressing the unquenched quantities in terms of the quenched ones. 
Subsections 3.2 and 3.3 then give the final answers for the formulas and pictures of the 
spectral densities at weak and strong non-Hermiticity, respectively, including the explicit 
examples of two degenerate and non-degenerate fiavours. Section 4 is devoted to a detailed 
analysis of the sign problem. In subsection 4.1 we first analyse the boundary, height, and 
frequency of the highly oscillating region of the spectral density. In subsection 4.2 we then 
define an observable which measures the severity of the sign problem and compute this 
observable for low and high density. We conclude in section 5. The derivations of several 
technical details are given in three appendices: Appendix A deals with the weak limit 
of the quenched spectral density of real eigenvalues, in appendix B we explicitly verify 
the decoupling of heavy flavours, and appendix C expresses the (partially) sign-quenched 
partition functions in terms of known spectral densities. 

2 The partition function and two limiting theories 

In this section we briefly introduce the chiral RMT and discuss its two different regimes of 
applicability in the physics of two-colour QCD. 

The RMT partition function at flnite matrix size and fixed topology v is given by^ 

-^Tr(PP^ + QQ^) 

Nf / \ 

xTTdet + (2-1) 

Here, the rectangular matrices P and Q are each of size N x (N + u) with real elements, 
without further symmetry restriction. The prefactor is chosen to ensure that4^^=°'^) = l. 
The parameters and mj denote the chemical potential and quark mass, respectively, 
for each fiavour / in the two-matrix model. The mapping of these RMT parameters to 
physical parameters is given in eqs. (2.4) and (2.6) below. The RMT in eq. (2.1) was 
introduced in [30] as a model for two-colour QCD at nonzero chemical potential. At = 
the matrix Q decouples, and the remaining one- matrix model equals the known RMT for 
two-colour QCD [41]. In eq. (2.1) the Dirac matrix is chosen to be Hermitian in the limit 
/xj — )• 0, hence we must take the to be purely imaginary in applications to lattice QCD. 
In [34] the above RMT was solved in the quenched case (Nf = 0) for finite N and in two 
different large- A'^ limits, i.e., weak and strong non-Hermiticity as explained in section 1. In 
[33] the precise limits in which the above RMT matches two-colour QCD with ^ 7^ were 
identified. 

^Note that our convention for the Gaussian weight is the same as in [34, 37] but different from [33]. 



zif^'"^(A.;{m}) = — I dP I dQ exp 
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In the first large- limit, the microscopic limit at weak non-Hermiticity (w), all chem- 
ical potentials are set equal (/ij = n) and are taken to zero when N ^ oo such that the 
product Nfj? remains finite. Likewise, the chiral limit is taken for all masses such that 
the product ^/Nrnj remains finite. By integrating out the random matrices, followed by a 
Hubbard-Stratonovich transformation and a saddle-point approximation, it was shown in 
[33] that the RMT partition function (2.1) maps to the static part of the chiral Lagrangian, 

Zi^^'"^(;u;{m})~ j dee"''ldUexp[^Nfi^TT(uiU^B^{UIU^)^B) 

U(l) SU(2Afj) 

+ \/iVReTr(e^^/^/?7/C/^iT"f)j , (2.2) 

which provides the dominant contribution in the e-regime of chiral perturbation theory. 
Here, we have introduced the 2Nf x 2Nf matrices 

-("xi). -(rr)' 

where t = diag(mi, . . . ,m7Vj), B is the baryon charge matrix, and T is the mass matrix. 
Comparing eq. (2.2) with the static limit of the chiral Lagrangian [9] we obtain the following 
mapping of RMT quantities to physical quantities, 

A' = 2iV/i2 = 4f4^y,F^V^ , (2.4a) 
m/ = 2\/]Vm/ = 2m/_physG'V4 (2.4b) 

with the physical parameters fiptys and mj^phys, the four- volume V4, and the low-energy 
constant G (in the notation of [9]). In particular, for equal masses these can be related to 
the pion mass ttt-tt and the low-energy constant F hy m = 2m'^F'^V4, where we have used the 
Gell-Mann-Oakes-Renner relation. The scaling of eq. (2.4b) at small chemical potential 
does not come as a surprise since it also holds true for vanishing chemical potential [42]. 

In the second large- limit, called strong non-Hermiticity (s) in [34], a map to a 
completely different regime at high density was found in [33]. In this limit the RMT 
chemical potential is not scaled, i.e., = 0(1), and it was found for an even number Nf 
of flavours that 



zi^-^'''-°\^ = l;{m}) ~ / / dUdV e^p -^ReTr{e'^'hUIU^t'^V* IV^) 



(2.5) 



U(l) S\J{Nf) 

in the sector of zero topology. Here we have to identify 



^) = ^) = :^^4AVj phys , (2.6) 

where A is the BCS gap. Note that we have rhf = nif, i.e., no rescaling by a power of N 
is necessary in this regime [34]. 

Whilst in [33] this mapping was found at maximal non-Hermiticity {fi = 1), we will 
show below that the same result holds for < // < 1, both for the partition function and 
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for the eigenvalue density correlation functions, subject to a trivial rescaling. We also note 
here that the two different limits in ^ can be formally related. By taking the limit /i — )• oo 
of the results in the weak large- iV limit (and suitably rescaling the eigenvalues and masses), 
the results in the strong limit can be recovered. 

The group integrals in eqs. (2.2) and (2.5) are in general very difficult to evaluate 
explicitly. However, in [37] a Pfaffian expression was computed for the RMT partition 
function (2.1) for finite N and any number of flavours. Combined with the large- limit 
of the building blocks of the Pfaffian expression, the kernel from [34], and the limiting 
individual skew-orthogonal polynomials of even degree, we find 



-)[Nf/2] 



Aiv,({m2}) 



Pf 

l<f,9<Nf 



Pf 

i<f,g<Nf 







2 



Nf even . 



Nf odd (w) , 



(2.7) 



where [x] denotes the integer part of x and the Vandermonde determinant is defined as 
'^n{{z}) = Ylkyii^k — zi). Note that the limiting partition function for the strong case 
only exists when there is an even number of flavours. The building blocks in the weak limit 
{w) are the quenched limiting kernel 



1 



2567r/i 



r,2 



(2.8 

and the limit of the quenched even skew-orthogonal polynomial g^v given in eq. (A. 4), 



lim 



1 



Z \^/2 



JV^^io (4iV)'^/2iV! (4iv) 

-'^'/V.(V^). 



(In 



In the strong limit (s) we have^ 



/C,(n, v) = '^{u-v) e-''-("+-) (2r/+V^) , 



where we have defined 



1 



(2.9) 



(2.10) 



(2.11) 



These objects will also appear as building blocks in expressions for the densities of Dirac 
eigenvalues. They are related to the quenched kernel K.^ (given explicitly in eq. (A. 3)) as 

u/2 

fC^{u,v) = lim vT^TTTT ( , .t-. ^ 1 fC 



AT^oo (4iV)2 V(4A^)2 
lCg{u,v) = {uvY^"^ lim ]C]y{u,v 



N 



/ n V \ 

viiv' m) 



/2NtJ. fixed 



(2.12) 
(2.13) 



*In order to emphasise the similarities between the weak and strong case, we have incorporated (by 
convention) an extra factor of {uv)"^"^ into the hmiting strong kernel K.^ compared with [34]. Corresponding 
factors have been removed from the limiting weight functions in eqs. (3.31) and (3.36) below. All observables 
remain unaffected by this choice. 
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We also note that the building blocks are directly related to the one- and two-flavour 
partition functions, 



'^(^;mi,m2)- ' , (2.14) 



2^ 

mj — mo 

Z):,'^-"''{(i-m)^q^{m^). (2.15) 



From the explicit knowledge of the partition functions above one could in principle 
derive detailed sum rules for the Dirac operator eigenvalues, as was suggested in [13, 43]. 
However, we will be able to give more detailed unquenched density correlation functions 
in the following section, from which such sum rules also follow. 

3 Dirac eigenvalue correlation functions 
3.1 Results for finite AT 

In the following we will be interested in finding the correlations of the eigenvalues of the 
Dirac operator. Within RMT this operator is given by the matrix 



V{,)^[ + 0^ ^ ), (3.1) 



with the Gaussian weights specified in eq. (2.1). In our convention, the Dirac operator is 
symmetric at // = 0, with the eigenvalues lying on the real axis. Since the Dirac operator in 
Euclidean field theory is actually anti-symmetric, it is necessary to multiply all eigenvalues 
by i when mapping from RMT to QCD. 

Finding the spectrum is equivalent to finding the eigenvalues of the Wishart matrix 
AB^ , as can be seen from the characteristic equation 

N 

= det[Al27v+i. - = A" det[A2ljv - AB^] = A" JJ(A2 - A^) . (3.2) 

i=i 

Because the matrix elements of A and B are real the solutions of the characteristic equation 
in A^ are either real or come in complex conjugate pairs. Consequently the determinant of 
the Dirac operator is always real, but not necessarily positive. Hence, for two-colour QCD 
there is no phase problem but a true sign problem. 

Switching back from Wishart to Dirac eigenvalues Aj we therefore have three possibil- 
ities: The Aj can be real (A^ > 0), imaginary (A| < 0), or come in complex quadruplets 
(±Aj,zbA*). In addition there are generically eigenvalues equal to zero. 

For simplicity we will assume that iV = 2n is even. Because we will subsequently take 
the large- A'^ limit this is not important.^ The partition function (2.1) can be expressed in 



^Results for odd A'^ are also known, and it is expected that they will lead to the same large-A'^ spectral 
densities. 
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terms of the Wishart eigenvalues Zj = A^, and the most compact form derived in [37] reads 



\Zi,...,ZN) 



Z^^^^-^^\^;{ra})^\{jd'z,p\,''^\^ 

Nf N N/2 

~ H"^/!! / d^^u\{F^''^\z2j-i,Z2j) ^n{{z]) , (3.3) 
/=i fc=i I i=i 

where in the second hne we have ordered the eigenvalues. This is achieved by the anti- 
symmetric bi-variate weight function defined as 

Nf 

F^^f\zi,Z2) = llizi- m}){z2 - mj){2ig{zi,Z2)[Q(lmz,) - Q(lmz2)] 5\z2 - zl) 
/=i 

+ /i(zi)/i(z2)'5(Imzi)(5(Imz2) sgn(Re Z2 — Re . (3.4) 

In contrast to [37] we have made the mass terms explicit in F^^f \ It contains the quenched 
weight functions 

Mz)^2|z|^/V-^i^|(r?+|z|), (3.5) 
g{zi,Z2) = 2|ziZ2|"/2g''-(^^+^^) p'^^-^lti4+-4)~V^tj^^_ {2^ltziZ2) erfc (r^+^t\z2 - zi|) , 

(3.6) 

where K and erfc denote the modified Bessel function of the second kind and the com- 
plementary error function, respectively. We note that in the limit of real arguments x, 
\im.z^.2^xg{zi,Z2) = h{x)'^. 

We are now able to give the result for all density correlation functions. The fc-point 
density correlation function for Wishart eigenvalues is defined in the standard way. For 
finite (and even) N, it reads 

Rl ^'{zi,...,Zk) 



Pf 

l<i,j<k 



d^zlc'^^^\z„z)F(''f\z,z,) 



I <fz,^,-.. j^<fzN P\^^\{z}) 



d'zfC^^f\z.,z)F^'^f\z,z,) 



fW{z,,z,)- I d^zF^^f\z^,z) 

c 

X / d^z')C^^''\z,z')F'^''f\z',z,) 



■ (3.7) 



It is given by a Pfaffian of the ordinary, 2k x 2k matrix composed of the matrices inside 
the square brackets, containing the unquenched kernel k!"^^\u,v) and single or double 
integrals folded with the anti-symmetric weight function F^^f \u,v) from eq. (3.4). Thus 
for fixed k there will be various contributions to r\^^\ beginning with one having only k 
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complex eigenvalues, then ranging over all mixed possibilities, to one having only k real 
eigenvalues. 

The unquenched kernel can be expressed as follows [37] in terms of the quenched kernel 
ICn (given explicitly in eq. (A. 3)) at finite (and even) N, 



Pf 



f^N+Nfi^^^) f^N+Nfiu,ml) 



7v+7v/"^/,^i) lCj^^j^^{mj,z) /Cjv+jv^Xm-p^is). 



WjU{u-m?j){z-m)) Pf /C^_^^^(my 



{31 



This expression is valid for an even number Nf of flavours, and even N. A similar expression 
holds for odd Nf as given in eq. (5.32) in [37] and also includes skew-orthogonal polynomials 
of even degree, as given in eq. (A. 4) (cf. eqs. (3.25) and (3.26) for Nf = 1 as an example). 
Because of eq. (2.7) the unquenched kernel can also be written as the ratio of two partition 
functions, for both even and odd Nf. 

Let us also comment on the two different large- limits that can be obtained from 
eqs. (3.7) and (3.8). In the strong limit for an even number of flavours, all fc-point corre- 
lation functions follow unambiguously by replacing the kernels in eq. (3.8) with the one in 
the strong limit, eq. (2.10). At weak non-Hermiticity the large- A^ limit is more involved, 
as was already noted in [34], due to the non-interchangeability of the limit and the inte- 
gration.^ Considering only complex eigenvalues this is not an issue as all integrations are 
"killed" by the delta functions in F^^f^ from eq. (3.4), and the A;-point densities follow by 
using the kernel in the weak limit, eq. (2.8), in eq. (3.8). Once correlations of real Wishart 
eigenvalues are computed the non-interchangeability becomes significant, as is explained 
for the quenched building blocks in appendix A. For the spectral density of real eigenvalues 
we explicitly evaluate the limiting expressions for Nf = 0, 1, 2, and then a higher number of 
flavours readily follows. Higher fc-point correlation functions in the weak limit that contain 
real eigenvalues will also involve double integrals (see the lower right block in the Pfaffian 
(3.7)) which may require further analysis. 

Although the structure of the most general fe-point function with arbitrary Nf is clear 
from the results above, in the following we will focus mainly on the most useful examples, 
that is, the spectral densities of complex, real, and imaginary Dirac eigenvalues with few (or 
no) flavours. They can be obtained from the spectral density (i.e., the one-point function) 
of the Wishart eigenvalues, which follows from a single kernel and is given by 

<^)(zi)= / d^zX:!f^)(zi,z)F(^/)(z,zO = M''^''^^(.i) + 5(yi)M''^«(xi), (3.9) 
JC 

where we denoted zi = xi + iyi. The delta functions in the first and second line of F^^f^ in 
eq. (3.4) lead to the split into separate densities of complex and real Wishart eigenvalues. 



6 



This fact has nothing to do with the sign problem as it also occurs at ^ = [41] . 
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The densities of complex, real, and imaginary Dirac eigenvalues are then obtained from 



(Nf,l 
Pi 




= 4| 


z\ 




(Nf,t 
Pl 




= 2| 


X 




(Nf,iK 
Pl 


\iy) 


= 2| 


y\ 





(3.10) 

To clarify the difference to the non-Hermitian RMT for three-colour QCD with /i 7^ 
[26, 27] let us briefly summarise the distinctive features of the spectral density for complex 
Dirac eigenvalues in the present case: 

(A) p^^''^\z) has no imaginary part. 

(B) p[--^\z) = pf-^\z*) = pf-^\-z*). 

(C) p^^''^\z) becomes zero for real and imaginary z. 

(A) and (B) follow from the reality and chiral symmetry of the Dirac determinant, while 
(C) is a result of the eigenvalue repulsion: If z approaches the real axis, z and z* come close 
to each other, but the corresponding probability is highly suppressed by the Vandermonde 
determinant An{{z'^}) of Dirac eigenvalues, and the same statement holds for z and —z*. 
(A)-(C) hold for any finite A^, and as we will see also in the large- A^ limit. 

We now turn directly to the results in the two different large- A^ limits. For some details 
in the quenched case, see also [37]. 

3.2 Eigenvalue densities at low density 

In this subsection we give the results valid for the microscopic limit at weak non-Hermiticity 
{w). In this limit we rescale the chemical potential, the masses, and the Dirac eigenvalues 
in the following way (see eq. (2.4)), 



p' 


= 2Np'^ , 




= 2\/]Vm/ 







(3.11) 

Under these scalings, and at large A^, the microscopic densities for Wishart (squared) 
eigenvalues are given by 

We can then switch back to Dirac eigenvalues using the mappings in eq. (3.10). Because 
R^^^'^\—x'^) 7^ R^^^'^\x'^), i.e., the function is not symmetric in general, the densities of 
real and imaginary Dirac eigenvalues will also differ in general: p[!^''''^\s,) 7^ pl!^^'^^\iO- 
Later in the strongly non-Hermitian limit a symmetry between the real and imaginary 
densities will emerge. 
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We will first give the quenched quantities. As was shown in [34] from eqs. (3.9) and 
(3.4) the density of complex eigenvalues is given simply by the rescaled weight function 
ffio(C*^)^^) multiplied by the weak limiting kernel in eq. (2.8), 

p(f/=°'^)(e) = m\'uc',e)icAe,c') , (3.13) 

where 

5«,(^,^*) = -5^(^*,^) = 2zsgn(Imz)^lim^^ J 5 (^|^, (3-14) 

= 4« sgnlim z e*'^ / — e ^^f" ^^K;±\ j\z\ eric — ^ imz 

^ ^ ^ t 2 V32/i4i I ; V4/i2' ' 

The weak limit of the density of real and imaginary eigenvalues is not quite so straightfor- 
ward. From eqs. (3.9) and (3.4) the density of Wishart eigenvalues e!\^^'^\x) is given by 
a real integral over the finite- kernel. It turns out that it is not possible to exchange the 
limit N ^ oo with this integration, a fact that already occurs for = and Nf = in 
[41]. In appendix A we derive the new quenched result 

pi^f=''(^^^\^) = -2\^\G^(^e,e), (3.15) 

where 

G'-C^'^O ^ - r ^itln ii^-'^'' f dy+ i /^i./2 r ^^y] ^^i^^v)hM (3.16) 
[sgn{x')Yi^ \ 7_oo [sgn(x')]''/2 Jo 



I 

327^ 



Here we have defined the rescaled real weight function 

u/2 

nA^^oo4nWx\J '"Kan J ^ ""2 \^8/i' 



- J- ^ (S^y '= (^) = ( M I (3.17) 



and used the exponential integral, which is defined as 

poo 

En{x) = / dtt-''e-''K (3.18) 



Equation (3.15) is plotted in figures 1 and 2 for jl = ^J{).2 and 1, each at = and 2. We 
observe that the presence of exact zero modes leads to a depletion of the spectral densities 
near zero. 

We now turn to the unquenched densities, where we will give two examples which 
are both new: two generic flavours {Nf = 2), including the degenerate mass limit, and a 
single flavour {Nf = 1). We start with a general remark: when expanding the Pfaffian in 
the numerator of the unquenched kernel in eq. (3.8), we find that it always contains the 
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Figure 1. Quenched spectral density of real (blue solid line) and imaginary (red dashed line) Dirac 
eigenvalues on the positive half-line in the microscopic limit at weak non-Hermiticity for /t = •\/0.2 
at = (left) and = 2 (right). 



iNf=0, (i) 




{Nf=0, {i 
Pw 




Figure 2. Same as figure 1 but for /i = 1 at = (left) and v = 2 (right). 



quenched kernel )CN+Nf{u,v) times a Pfafiian Pf [/Cat+at^ (m^, m^)] that is the same as the 
Pfaffian in the denominator and hence cancels; the remaining product in the denominator 
is cancelled by factors in the weight function. Therefore, on inserting this into eq. (3.9), 
the unquenched density is seen always to be of the form of the quenched density plus a 
correction term, 

(3.19) 



i?i^^=°)(z) + Ai?S^^)(z) 



for both the complex and the real Wishart eigenvalues. This argument holds for both an 
even and an odd number of flavours. 

The two-flavour density of complex eigenvalues is therefore obtained in a straightfor- 
ward manner, 



(3.20) 



with the quenched density from eq. (3.13) and the limiting weight function (3.14). It is of 
the same form as for the symplectic ensemble relevant for adjoint QCD [29]. In the limit 
of degenerate masses we simply have to Taylor-expand the ratio of quenched kernels in the 
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second line, where the denominator vanishes hke rhi — m2. For convenience we also give 
the limiting expression, 



(3.21) 



where 



+ Ji,(sVa;)(z^J,.+i(s^) -s^J,,(s^)) . (3.22) 



Looking at eq. (2.14) and noting that is the derivative of /C«, with respect to the second 
argument, eq. (3.22) with x = y also provides the mass-degenerate two-flavour partition 
function, 

^{Nf ^.i')^^. ^ / f/s53g-2/is ^^J^^-^(^5^) J^_-|^(si72,) _ J^(5l7i)2^ ^ (3.23) 



after using an identity for Bessel functions. The two-flavour density of the real and imag- 
inary eigenvalues can readily be obtained from the above building blocks, using eqs. (3.8) 
and (3.9), 

(Nf=2,m), . _ m,{^)M)..^ , ^1 ^1 /C^(g^ mi)G«,(mf , -e^) - /C^(g^ mf )G^(m^, g^) 

Pw {z,, rni, m2) — [^q,) i" ^|?| 



(3.24) 

with Gui deflned in eq. (3.16). However, there is no simple expression for the degenerate- 
mass limit in this case. 

Our second unquenched example is the case of a single flavour {Nj = 1). The density 
of complex eigenvalues reads 

(3.25) 

which includes the second building block eq. (2.9), 



For the density of real and imaginary eigenvalues we again encounter the problem of non- 
commutativity of the large- limit and integration. By handling this issue similarly to the 
quenched case (see above), we obtain 

pT^=''^'>^\i;rh) = pL^^=°'(-^^)(g) + 21^1 ^-^^')^-^^''^') 7 ^ (3.26) 
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Figure 3. Quenched (left) and Nf = 1 (right) spectral density of complex Dirac eigenvalues at 
weak non-Hcrmiticity, both for fi = \/0.2 at = 0. The right figure is in the chiral limit (m — 0). 




Figure 4. The Nf = 1 spectral density of complex Dirac eigenvalues at weak non-Hermiticity for 
/t = 1.8 (left) and jj, — 2.5 (right) at z/ = in the chiral limit (m = 0). Here and in the following 
the peaks of the oscillations are clipped. 



where 



hw{x) 



n /-x 



[sgu{x)Yi^ y V J-oo [sgn{x)YI^ Jo J 

The proof follows along similar lines to that of G^(m, u). 

In figures 3-5 we show plots of p^w^ from eq. (3.13) and p"w^ ^'^''(C;"^) from 

eq. (3.25) for several values of ji and rh. (Due to the symmetry p^w^''^\C) = Pw^''^\0 = 
Pui^^'^\—£,) we only present the result in the first quadrant. The quark masses are taken 
to be imaginary as our Dirac operator is symmetric at p = 0.) For small fi (figure 3), 
the effect of the massless flavour is to deplete the spectral density near the origin. As fi 
increases (figures 4 and 5), qualitatively different effects emerge. In the region between the 
location of the mass (in this case, the origin) and the edge of the spectrum, there appears 
a domain of strong oscillations with the shape of an ellipse that grows in its peak size and 
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Re[e] Re[^] 

Figure 5. Left: Nf — 1 spectral density of complex Dirac eigenvalues at weak non-Hermiticity 
for jj, = 6 and rh = 20i at = 0. Right: Envelope of the oscillating region given by eq. (4.14), see 
section 4.1 for more details. 




Figure 6. The Nf = 2 spectral density of complex Dirac eigenvalues at weak non-Hermiticity for 
fl — 6 and v = 0, for mi = m2 = lOi (left) and rhi = lOi and m2 = 70i (right). In the right figure 
only the positive part of pin ^~ ' is shown for better readability. 

in the number of ripples as fi gets larger (see section 4.1 for a detailed analysis of these 
oscillations). 

This change in the spectrum signals the emergence of a severe sign problem. For Nf = 1 
the Dirac determinant is real but not necessarily positive, hence the spectral "density" 
could take negative values. The structure of the oscillating domain, to be analysed in more 
detail in section 4.1, is quite similar to the corresponding one in the non-Hermitian RMT 
for three-colour QCD with /i 7^ [27], despite the difference in symmetries between two 
and three colours. 

Another important remark is in order. When rh becomes so large that it is located 
outside the support of the quenched spectrum, the spectrum shows no change as compared 
to the quenched case. This is consistent with the physical expectation that a heavy flavour 



- 15 - 




Figure 7. The Nf = 1 spectral density of real (left) and imaginary (right) Dirac eigenvalues at 
weak non-Hermiticity for /t = 3 and i' — 0. Left: rh — O.li (black), m ~ 2i (blue), rh — 3i (red), 
and quenched case (dashed). Right: m = 12i (black), m = 15i (blue), m — 2Qi (red), and quenched 
case (dashed). Note that the zeros coincide with the mass values. 



should effectively decouple from the rest of the system. More importantly, it suggests that 
the sign problem will be milder when the quark mass is larger than a certain critical value, 
set by the width of the quenched support. An essentially identical conclusion was reached 
for the three-colour RMT (see, e.g., [44] for a review), again despite the difference from 
our two-colour model. 

Next, we turn to the plot of pw ^ ' mi, 7712) from eq. (3.20) for two sets of masses 
in figure 6. For degenerate masses (left), there is no sign problem and no change in 
the spectrum except for a depletion of eigenvalues near the location of the masses. This 
situation is qualitatively very similar to the symmetry class for adjoint QCD (which has 
no real or imaginary eigenvalues), see the figures in [28]. For widely different masses, the 
density develops a domain of oscillations, which resembles the A'^^ = 1 case. This implies 
the decoupling of one flavour to become heavy. In section 3.3 it will be shown that a similar 
oscillating spectrum emerges at high density for well-separated masses for Nf = 2. This 
was to be anticipated from the connection between the weak and strong large- limits (see 
the comment in the paragraph following eq. (2.6)). 

In figure 7 we show plots of pH^^ ^' ^''^^^ (^; m) from eq. (3.26). The effect of a dynamical 
flavour is always stronger for the imaginary eigenvalues than for the real eigenvalues. The 
reason is that the partition function and densities vanish when the imaginary eigenvalues 
equal the (in our case) imaginary masses, see eqs. (3.4) and (3.3). Thus the imaginary 
density crosses the horizontal axis at the value of the mass.^ Interestingly the dynami- 
cal flavour brings about an oscillatory behaviour in the real density. For both real and 
imaginary densities, one can clearly observe a convergence to the quenched density as the 
mass is increased, which leads to an effective quenching (see appendix B for an analytic 
confirmation of this fact). 

^Using eq. (C.7) in the appendix and the positivity of the partition function for imaginary masses, one 
can prove for general Nf that the imaginary density cannot cross the horizontal axis at values away from 
the masses. This explains why a rapid oscillation as in figure 5 does not occur in the imaginary density. 
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Figure 8. The Nf = 2 spectral density of real (left) and imaginary (right) Dirac eigenvalues at 
weak non-Hermiticity for /i = 3 and ly — 0. Left: (7711,7712) = {0.1i,0.1i) (black), (O.lijlOi) (blue), 
(0.1«,307) (red), and Nf — 1 case with ttt, = Q.li (dashed). Right: (7711,7772) = {12i,12i) (black), 
(127,207) (blue), (127, 30«) (red), and Nf — 1 case with 777 = 127 (dashed). Note that the zeros 
coincide with the mass values. 

In figure 8 we show plots of plf^^ mi, 7712) from eq. (3.24). We illustrate the 

effect of dynamical flavours by changing one of the masses while keeping the other one 
fixed. For both real and imaginary densities, one can clearly observe a convergence to the 
Nf = 1 density as one of the masses is increased. Note also that the spectral density for 
degenerate masses is positive, as it should be, because there is no sign problem there. 

Because of the issue of non-interchangeability of integration and large- limit in the 
weak limit, as well as in order to better understand explicitly the relation between theories 
with different flavour content, we have checked in appendix B that sending a subset of 
masses to infinity leads to the decoupling of these. 

3.3 Eigenvalue densities at high density 

In this subsection we will give the results for the spectral density in the large- limit at 
strong non-Hermiticity, which corresponds physically to the high-density regime. In this 
large- limit the rescaling of the parameters is different from the one in eq. (3.11) and is 
given by (see also eq. (2.6)) 



mf = mf , 
e = A, 



(3.28) 



i.e., the masses and the Dirac eigenvalues, as well as n, need no rescaling at all. For the 
microscopic densities of Wishart eigenvalues we have 



i?f^'^)(z) ^ lim <^'^)(.) 

A''—!' 00 



lim <^«(x) 



(3.29) 



and the densities of Dirac eigenvalues are again obtained using the mappings in eq. (3.10). 
We again begin by quoting the quenched densities from [34], which will constitute the 
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building blocks for the new unquenched results. We have 

^=°'^^(e) = m\'9sic',e)ics{e,c') , (3.30) 

where the complex weight in the strong limit is defined as 
9s{z,z*) = -gs{z*,z) = 2i sgn(Imz) j^g{z*,z) (3.31) 

/•CO Jj. 

= 4isgn(Imz)e2^-^^^ / -e-''+*(^'+^*')-si^. {2r]lt\z\'^) edc {2ij+Vt\lm z\) 
Jo t ^ 

with g from eq. (3.6). It can be seen from that equation that the cases of general ^ < 1 
and of maximal non-Hermiticity with = 1 (and r/+ = | ) are related simply by a rescaling 
of the complex eigenvalue — )• The same is true for the densities of the real and 

imaginary eigenvalues, 

^{JV,=o,WM)(^)^2|e|^2K.(,?+|e|2)(^"rfx'|e'-x'| 2i^.(r,+ |x'|)/.(2r?+e^^) (3.32) 

+ r dx'\e-x'\ 2K|(r?+|x'|) J,(2r?+ex^)) . 

Unlike in the weak limit, in the strong limit no complications arise from an interchange of 
the large- limit and the integration. 

The example of two non-degenerate fiavours is now straightforward, and we only quote 
the result which includes eq. (3.30) as the first term, 

{Nf=2,C), . . . 
Ps (4;mi,m2) 

= Pi»«(f) + 4KP^.(r^{^)^^«'^"^'^^«"^5'-'^f'■"'''^'«""^^' . (3.33) 



It is identical to the form of the weak result in eq. (3.20) apart from the kernel in the strong 
limit from eq. (2.10), 

3 

It is easy to obtain the limiting result for two degenerate masses which we also give here 
for completeness, 

Ps i^; mi, 1712) 

= ^=°'"^(0 +4|el^^.(r^af^^^^^^ 

8TTl^{2r]^mf) L 



+ (r - mi){C' - mf)^ /,(27?+mir)a.+i(2r/+miO " (^ ^ D • (3-34) 

nil ^ ' 

Similarly, the extension of eq. (3.32) to A'^^ = 2 reads 

pf^=^'«^)(e;m,,m2) = ^=°'«^^(0 + 2|e|/^s(e^) / dx /..(x) sgn(e^ - x) 

Jr 



/C,(C^ m^)/C,(x, mj) - /Cs(g^ mf )/Cg(x, m|) 
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Figure 9. Quenched (left) and Nf — 2 (right) spectral density of complex Dirac eigenvalues at 
maximal non-Hermiticity (/i — 1) for v — 0, with equal rhi — 1112 — 2i in the right figure. Note 
again the similarity with adjoint QCD [28]. 





Figure 10. The Nf = 2 spectral density of complex Dirac eigenvalues at maximal non-Hcrmiticity 
(fi = 1) for v = with (toi, 7712) = {2i, 5i) (left) and (toi, 1712) = (2i, 8i) (right). In the right 
figure the negative part of ps is cut for better readability. 



where the strong limit of the real weight is defined as 



hsiz) 



1 



h{z) 



\u/2 



i-z) 



u/2 



(3.36) 



The density ' 'i^; rhi, 1112) is plotted in figures 9 and 10. In the former the quenched 

result is also included for comparison. (Again, we only present the result in the first 
quadrant because of the symmetry pi^^''^\£,) = pi^^''^\^*) = Ps^^''^\—0^ we take 
the quark masses to be imaginary as our Dirac operator is symmetric at /i = 0. Since 
the difference between < p < 1 and p = 1 is trivial we set p = 1 in the figures.) We 
observe that the effect of dynamical flavours with equal masses is to simply generate a 
dip in the quenched spectrum. The effect of unequal masses is qualitatively different. For 
a sufficiently large mass difference, the spectrum exhibits a domain of strong oscillations, 
taking both positive and negative values, which signals the appearance of the sign problem. 
We caution that sending one of the masses to infinity does not yield pi^^ ^''^\^;m) but 
merely makes the oscillating domain extend to the entire spectrum.^ 



^This is not in contradiction with the usual decoupling of heavy flavours, see appendix B, since in the 
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Figure 11. The Nf — 2 spectral density of Dirac eigenvalues at maximum non-Hcrmiticity (p = 1) 
for (mi, 7712) — (2«, 3i) at = (left) and v ~ 2 (right). The densities of both real (blue full line) 
and imaginary (red dashed line) eigenvalues are displayed in the same plot for comparison. The 
latter again has zeros at the locations of the masses. 



The density ps ^ ' * ((*)■?; "^ii "^2) is plotted in figure 11. As in the case of weak 
non-Hermiticity, the density of imaginary eigenvalues goes through zero at = irhi and 
^ = 11712. However, unlike in the case of weak non-Hermiticity, no oscillation of the density 
of real eigenvalues is observed. A more detailed analysis of the sign problem will be 
presented in section 4. 

We now consider the case of four quark flavours, which is interesting since at high 
density Nf = 2 and Nf > 4 (even) have distinct chiral symmetry breaking patterns and 
a new class of Nambu-Goldstone bosons emerges for Nf > 4 [13]. With non-degenerate 
masses, we have to expand the PfafRan of a 6 x 6 matrix in order to write down the complex 
density explicitly. The result is 

Ps (C; {m}4) = ps (0 (3.37) 

ICsirhf, m|)/Cs(m|, m|) - Ks{m\,rh\)Ks{ml, rh\) + /Cs(mf , ?fi|)/Cs(m|, m^) 
X |/Cs(mi,^*^)/Cs(m2,C^)/Cs(m|,m4) - /Cs(mi, ^^)/Cs(m2, ^^)^s("^i, m|) 
- 1Cs{ml C^)lCs{ml ml)lCs{rnie) + /C.(m2, m2)/C,(m2, r2)^^(^2^ ^2) 

+ K.s{m\,C'^)lCs{rhl, ml)lCs{rh\, ^^) - /Cs(mf , m|)/Cs(m2> r^)^s("ii C^) 
+ K,{ml rhl)]Cs{ml r')/C,(mi ^2) _ /^.(rn?, e%s{rnlrnl)lCs{rnl C') 

+ )Cs{ml ml))Cs{mie)}Cs{ml C^) - /^(m?, m2)/C.(mi, e))Cs{ml f')} . 

Plots of pi^'^ ^''^\^; {rh}i) are shown in figures 12 and 13. It is seen that Nf = 4 
flavours roughly amount to two pairs of Nf = 2 flavours, each of which shows either a dip 
or a domain of strong oscillations. 

strong limit we are restricted to even Nf. 
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Figure 12. The Nf = 4 spectral density of complex Dirac eigenvalues at maximal non-Hermiticity 
(/X = 1) for mi ~ 7712 ~ = "tn^ — 3i at v — (left) and v — Q (right). 




Figure 13. Same as figure 12 but for toi — TO2 = 4i, 7713 = TO4 = 8z (left) and toi = TO2 = 3i, 
ma = 5i, 7714 = llz (right), both at = 0. Note that the scales for Re[^] and Im[^] are different. In 
the right figure the negative part of ps is cut for better readability. 

4 The sign problem 

4.1 Analysis of oscillations of the Dirac spectrum 

If we choose non-identical masses the spectral densities show interesting behaviour. In this 
section, we will consider in particular the case of Nf = 2 with rhi 7^ m2, and the case of 
Nf = 1 in the weak non-Hermiticity limit. As can be seen from a direct inspection of the 
expression (2.1) for the partition function, the presence of non-degenerate masses can lead 
to "probabilities" which, although still real-valued, may be negative in certain regions. This 
is perhaps not unexpected. However, what may be surprising is that there are regions of 
the complex plane in which the eigenvalue density function seems to oscillate from positive 
to negative values quite dramatically. We present here a straightforward explanation. 

For simplicity, we first consider the strong limit which was discussed in section 3.3. In 
the quenched case, the complex Dirac eigenvalue density is essentially constant everywhere 
apart from close to the real and imaginary axes, where it vanishes. As we discussed earlier 
in eq. (3.19), the unquenched density can be written as a correction to the quenched case. 
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{-2i) Im [iCsif M)^s{e\ ^l)} , (4.1) 



Prom eq. (3.33), the correction in this case is 

where the weight function (jsiz* , z) was given in eq. (3.31). We next insert the exphcit form 
of the strong kernel from eq. (2.10), and, after noting that a number of factors cancel, we 
get 



exp[-r?„(e^+r')] 



X {-2i) lm[{e - m?)(r' - ml)h{2l^+im{)U{2l^+Cm2)] . (4.2) 

Here, r/+ and ?]_ are as defined in eq. (2.11). So far everything is exact. However, we 
now make two approximations. First, we approximate 5s('C*^;?^) by using the following 
argument. For the quenched case, away from the axes, the Dirac density is essentially flat, 
i.e., it is equal to some constant which we could determine exactly, but here we shall simply 
denote it by ps- Therefore, from eq. (3.30), we have 



m'ue.nK-sie^n^ps. 

On inserting the explicit form of the strong kernel, and rearranging, we have 



(4.3) 



4|C|^ 



Svr 



*2\ 



PS 



(e2-e2)/.(2r?+|eP 
and so, away from the axes, 

(iv,=2,c),^.^ ^ , _ (-2i)/95 Im{(^2 _^2)(^*2 _^2)j^(2r?+emi)/,(2r/+e"i2)} 



(4.4) 



(^;mi,rfi2) 



-2 



[m 



I - ml) I,{2r^+mim2){e - i*^) Iu{2r,+ \m 



(4.5) 



Our second approximation is to replace each modified Bessel function with its large- 
argument asymptotic limit, namely 



h{z) 



^/2 



■KZ 



(4.6) 



for — f < arg(z) < ^, i.e., in the right half-plane. After simplification, this gives (in the 
right half- plane, and sufficiently far away from the origin and from the axes) 



iXps {t,;mi,m2) 



{-2i)ps e" 



-2r7+(mim2 + |5|2) 



{m{-ml){e -i*^) 
{-'ii)ps 



exp 



2??+ ((mi + m2)x - rhim2) 



{ml-ml){e -i*^) 

X Im |(^^ — m\){^*'^ — ml) ^ cos(2ry+(mi — rh2)y) + i sin(27/+(mi — m2)y)^ | , (4.7) 
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where we have written x = Re^ and y = Im^. We can regard this function as being the 
product of an envelope (the exponential), some polynomials in ^ and and an oscillatory 
function of y. Let us consider each in turn. 

To understand the behaviour of the exponential part, we can determine the contours 
in the ^-plane where this part takes constant values. We can equivalently perform this 
analysis on the exponent. So we look for the solution of 

(mi + m2)x — — mirh2 = k , (4-8) 

where k is some constant. Writing = + y^, and rearranging, we have 



(^^S^) '-k^r\ (4.9) 

This is, of course, the equation of a circle of radius r, centred on the point ^ = (mi + 7712) /2. 
The exponential function equals unity when k = 0, i.e., when the circle passes through the 
points ^ = mi and ^ = 1712- Crudely speaking, this marks the boundary of the region where 
the correction is "significant" when compared with the quenched case (where the density 
equals ps)- Points inside this circle have A; > 0, so the circle marks the boundary of a peak 
(rather than a dip) in the complex plane. The maximal value of this peak is exp(2r/+/cniax); 
where 

^.ax=(^^^)', (4.10) 

SO the peak grows very rapidly indeed with the difference between rhi and m2. Of course, 
our approximations are not valid at the centre of the circle since this lies on the x-axis, 
and so this argument is merely indicative of the orders of magnitude that we might expect 
to observe in the density itself. 

The polynomials (^^ — mi)(^*^ — m^) in eq. (4.7) result in a repulsion of eigenvalues 
from the points ibmi and ±m2, similar to what we observed in the case of identical masses. 

Finally, if we combine the oscillatory part of eq. (4.7) with the {C'^ — (,*'^)~^ factor, then 
we get an oscillating function of y which is symmetric (even), and which has a constant 
"wavelength" 

A= t^"" ^ (4.11) 
ry+ (mi - m2) 

In other words, this constitutes a set of ridges and troughs parallel to the real axis. 

Therefore, what we observe in the total density is a circular "peak" of fairly dramatic 
oscillations, but these are almost completely suppressed elsewhere in the complex plane. 
There is also some repulsion from the location of the masses, although the exponential 
effect causing the peak is generally more significant. The number of oscillations visible will 
therefore be given by 

-^oscillations ~ ; (4-12) 

TT 

and we can say that no oscillations will be seen if, roughly speaking, 

Imi — m2|<w — • (4.13) 
V V+ 
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These same conclusions can also be obtained directly from the group integral eq. (2.5). On 
making an ansatz for the matrices U,V € U(4) and keeping only the degrees of freedom 
that are expected to be relevant, a mean- field analysis reveals the same equations for the 
boundary and height of the oscillations. 

For the weak non-Hermiticity limit, a similar analysis can be performed, although 
the situation here is slightly more complicated because we need to use not one, but two 
approximations for the kernel, depending on the average of the arguments. 

When we have two masses, both located inside the strip that contains the quenched 
eigenvalues (and both taken to be purely imaginary), then we find that the situation is 
similar to the strong case above; there is a circular region of oscillations, with the boundary 
of the oscillating region passing through the two masses. Indeed, one can map the weak 
non-Hermiticity case onto the case of maximal non-Hermiticity (with /x = 1) by making a 
simple scaling of masses and eigenvalues by 1/2/i, and then the earlier results eqs. (4.11), 
(4.12) and (4.13) for the strong case also hold true here. 

If we have a single mass located within the strip at ih = irj, then we get an elliptical 
region of oscillations, given in fact by the same equation as for the corresponding situation 
in three-colour QCD [45, eq. (41)] where it was derived from an ansatz. 



A comparison of this formula with the exact spectral density is given in figure 5. Conversely, 
a single mass located outside the strip has virtually no effect whatsoever on the density. 

The Nf = 2 case with one of the masses considerably far outside the strip (in fact, 
beyond a (finite) critical point for which we can solve) is both qualitatively and quantita- 
tively very similar to the Nj = 1 case; we refer to Appendix B for further details on the 
decoupling as one of the masses is taken progressively larger. However, there is a third 
distinct intermediate regime where one mass is inside the strip, and the other is outside 
the strip, but below the critical point. Here, we find that the region of oscillations lies 
part-way between a circle and an ellipse. 

It would be interesting to study the effects of the oscillations analysed in this section 
on the chiral condensate, along the lines of [39, 40]. This is left for future work. 

4.2 Average sign factor 
4.2.1 Introductory remarks 

In this section we study the severity of the sign problem by looking at the average sign of 
the Dirac determinant. The definition of such an average is not unique, and we will use 
the following guiding principles to define an average that behaves reasonably: 

• The average should be bounded by 1 from above. 

• Starting at 1 for // = (or degenerate masses) the average should be a decreasing 
function of /i (or the mass difference). 





(4.14) 



(^(2/i2-^))' (i(2/i2-r?)) 
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Let us illustrate the problem of finding a proper definition using the case of three- 
colour QCD (QCD3), which has been well studied in the RMT framework [46-48]. If we 
denote the complex phase of det[P + m] by e*^, the authors proposed to compute 

/ 2i9\QCD3 ^ / det[V{fi) + m] \^^^^ 

^ " \det[P(M) + m]t/^^ ' ^^■'^> 

which gives the expectation value of twice the phase in the case that the determinant is 
complex. Had we chosen the opposite ratio on the right-hand side instead, this would 
naively have led to {e~^'^^)Nf However, for Nj = 1 we have in three-colour QCD 

/det[V + m]t _ (det[P + m]t)Q^3 



\ det[P + m] (det[P + m])Q^S ' ^^''^^ 

because the expectation value in denominator and numerator (also called characteristic 
polynomial) is equal and //-independent [49]. Obviously this would contradict the second 
of our criteria above. Note that in contrast using the definition (4.15) the ratio is not equal 
to 1, even for Nf = 1. 

Our definition of the average sign for two-colour QCD, which we will use throughout 
this section, is 

Nf 

p^''f\f,; {m}) = ( sgn n det[Z)(/i) + m/]\ (4.17) 
f=i II J^ii 

where the average is computed in a sign-quenched theory denoted by ||A^/|| (i.e., only the 
absolute values of the determinants appear in the measure). We recall that in the strong 
limit we always need Nf to be even. In the next two subsections we will consider two 
particular cases, the weak limit with one flavour, where 

pif^"^^(/i;m) = lim p^^f=^\fi;m) , (4.18) 

A''— >oo,woak 

and the strong limit with two flavours, where 

pij^f ^^(^; Tjj^^ 17^2) = lim p^^-f^'^\fi;mi,m2) ■ (4.19) 

A''— >cxD, strong 

We shall see explicitly from these two cases that our guiding principles are satisfied. 

We finish this subsection by commenting on possible bad choices for the average sign. 
Due to the obvious inequality 

( |det[P + m]| )^^^^ > { det[V + m])^^^^ (4.20) 

we have 

/ det[P + m]\ I \d(ii\D + m]| 



_ , > 1 ('4 21) 

\|det[P + m]|/^^,^i \ det[P + m] - ' 
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and hence a ratio of this kind, which would have been the natural generalisation of 
eq. (4.15), can be ruled out. A second choice that leads to an inconsistency is the follow- 
ing. Because the sign problem can be turned on by detuning the masses of two originally 
degenerate flavours one might consider 

/ det[V{fi) + mi] det[P(/i) + ma] , ^2) 



\ det[2?(;u) + mo 



2 



where in the denominator we have two degenerate masses mo • One caveat is that the choice 
of mo in the denominator is ambiguous. For small mi — m2 one might think that any of 
the choices mo = mi,m2, (mi + m2)/2, ^mima would be acceptable. However, using the 
latter choice we can show that in the strong limit this leads to a quantity that is always 
unity. Moreover, in the weak limit or for other choices of mo the quantity in eq. (4.22) 
goes to zero more rapidly for larger topology. This contradicts our expectation, which was 
confirmed in RMT studies of three-colour QCD [50], that a larger number of exactly zero 
eigenvalues should push the density away from the origin (where the oscillations are strong) 
and thus mitigate the average sign problem. 

4.2.2 Average sign at low density 

Let us consider the case Nf = 1 in the limit of weak non-Hermiticity. Although the sign- 
quenched partition function does not permit a field-theoretical interpretation like QCD 
with nonzero isospin chemical potential, it is still possible to compute it as a mathematical 
entity by RMT. In appendix C.3 we derive 

zi^^'="^"'^)(A;m) ^ ^^Inn ^^^(1 det[P(/.) + m]|>^^^„ 

^_^^,y,GJn^_ (4.23) 
/i^(m2) 

This result essentially follows from the quenched real spectral density. Combining it with 
eq. (2.15), 

„(Nf=l. v) , ^ ^ N / - 2\ 

Zw ' (/u;m) ~ Qwym ) , 
we obtain the average sign factor in the weak limit, 

= ~ "AeAV^O.(™^™=) ■ <^-^^' 

The normalisation of ^\fi;m) can be uniquely determined from the requirement 

pl^^='\0,m) = l. 

In figure 14 we depict how p[!^^ (Ai ™') depends on fi for several values of m and u. We 
observe that the onset of the sign problem is delayed as m or is increased, a feature which 
is in common with three-colour QCD [46, 50]. In figure 15 we plot the average sign against 



a rescaled variable \xl ^mj'l. In physical units, it is equal to /i/(mjr/2), where m^r denotes 
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Figure 14. Average sign for Nf = 1 at weak non-Herniiticity. Left: m ~ (blue line), m = 6i 
(red thick line) and rh — 15i (black dashed line), all for 1^ = 0. Right: u — (blue line), — 10 
(red thick line) and v = 20 (black dashed line), all for rh — 0. 



Pw 



(^f=^)fr..^^ 1.0 




2.0 7^ 



Figure 15. Average sign for Nf — 1 and = at weak non-Hermiticity. The curves correspond 
to TO = 4i (black), to = 8i (blue), m = I6i (red), to — 32i (magenta), and m = 64i (orange). 



the mass of the Nambu-Goldstone (NG) bosons in the vacuum.^ Interestingly, the sign 
problem looks almost absent for fi < while it deteriorates rapidly for fi > In 

particular, we observe a convergence of the curves to a step function in the thermodynamic 
limit. This apparent jump of the average sign is quite intriguing, compared to the average 
phase factor in the microscopic limit of three-colour QCD [46, 47]: the latter changes 
smoothly from 1 to in the limit m — )• oo. 

Since the /i-dependence of zi^-' ^' '^^ is explicit in eq. (2.15) (being of oscillatory nature, 
see eq. (2.9)), it must be zi,^^ M''^) ^^^^ responsible for the apparent jump of the average 
sign. While it seems that the Nf = ||1|| theory is similar to the Nj > 2 theory in that both 
undergo a transition at ^ = m^/2, see, e.g., [9], a comparison of the order of the transition 
needs further investigation. 

Finally we comment on a lattice simulation using one flavour of staggered fermions 
in the adjoint representation of SU(2) [7, 8]. This theory belongs to the same symmetry 
class as Nf = 1 two-colour QCD with fundamental fermions. In [7, 8] this theory was 
simulated by two different algorithms, Hybrid-Monte-Carlo (HMC) and Two-Step-Multi- 

^Note that no NG mode appears in Nf = 1 two-colour QCD. The above m-^ refers to the mass of the 
NG modes that appear in Nf > 2 two-colour QGD. 
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Boson (TSMB). The HMC algorithm is non-ergodic in this case and actually simulates a 
different theory, namely the one-flavour sign-quenched theory. The authors of [7, 8] report 
that the sign-quenched theory simulated by HMC seems to undergo a phase transition at 
/i = while the correct one-flavour theory simulated by TSMB exhibits no singularity 

at all. These results are consistent with ours. However, the observables considered here 
and in [7, 8] are different so that a direct comparison is difficult. 

4.2.3 Average sign at high density 

Let us now consider the case Nf = 2 in the limit of strong non-Hermiticity. This theory is 
free from the sign problem if and only if the masses are degenerate. Here we are interested 
in the severity of the sign problem for \rhi — nonzero. In appendix C.4 we derive 

2i^^""^"'''^(^;mi,m2) = Hm (| det[P(/i) + mi] det[P(/x) + ms]]) 

Af— >c«, strong \ ' / Nf=0 



fs{ml,ml 



— 7712 

where^*^ 

fs{rnl,ml) = ICs{ml,ml) sgn{ml - ml) 



(4.25) 



+ I dx j dx hs{x)hs{x') sgn(mi — j;) sgn(m2 — 2;') — sgn(7n2 — x) sgn(rni — x') 

K.s{rh\, x)ICs{iTi2, x) — -ICs{iTi1, rn^)K,s{x^ x) (4.26) 

with the kernel in the strong limit from eq. (2.10), 

3 

lCs{u,v) = |±(n-^;)e-''-("+^)/.(2r?+V^) . (4.27) 

These results are obtained from the two-point correlation function of real eigenvalues in 
the strong limit. Combining eq. (4.25) with eq. (2.14), 

„{Nf=2,v), „ , . ICsi'ihj,ml) 



mj — 1712 



we obtain 



Ps {fi;mi,m2)= (jv.-ll2ll = sgn - m^) ^ 2^ ■ (^-^S 

This expression satisfies the normalisation condition lim pi^^ 7111,1712) = 1. By 

m2— >rhi 

explicit calculation a scaling law is verified. 



p. 



^\pL-mi,m2)=pf' ^\l;^/2^mi,^/2^m2), (4.29) 



^"in eq. (4.26) the squared masses are assumed to be nonpositive since this is the case relevant for two- 
colour QCD. If one or both of the arguments of fs are positive we obtain a more complicated expression 



for fs{xi,X2) = limiv^oo(a;ia:2)''' fN{xi,X2) 
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Figure 16. Average sign factor for Nf = 2 at maximal non-Hermiticity (// — 1) with rfii — i (left) 
and fhi — Ai (right), at = (blue line), v = \{) (red thick line), and v — 2Q (black dashed line). 



hence we can assume = 1 without loss of generality. 

In figure 16 we plotpf^=^)(l ; mi, 777,2) as a function of 7772- The results indicate that 
larger topology mitigates the sign problem but becomes less effective for heavier masses. 
This phenomenon is easily understood if one recalls that a large v that causes a depletion of 
eigenvalues near the origin can smoothen the strongly oscillating domain of the spectrum if 
the domain is close to the origin; since the domain moves away from the origin for heavier 
masses, the smoothing by u becomes less effective. 

As a result, the average sign turns out to be essentially independent of topology for 
larger masses, where the sign problem becomes severe for \5rh\ > 2.5, as can also be seen 
from figure 16. In physical units we have, using eq. (2.6), 

|*,„„,„|>^. (4.30) 

We remark that the oscillations of the complex spectral density start if \5rh\ > ^/2^T ~ 2.5 
as seen from eq. (4.13) sX fi = 1. It is reasonable that these two values of |5777| are close to 
each other. 



5 Conclusions 

In this paper we have solved a random two-matrix model as a mathematical model for 
unquenched two-colour QCD with non-vanishing chemical potential. Our model has very 
interesting features. First of all it describes two very different physical situations, at low 
and at high density, in the approximation of a static chiral Lagrangian. Despite the sign 
problem the unquenched random two-matrix model can be solved exactly, and the Dirac 
operator has very distinct real, imaginary, and complex eigenvalue densities. The sign 
problem can be switched off or on at fixed chemical potential by tuning, say, the masses of 
two quark flavours to be degenerate or distinct. 

We have given many explicit examples for the various spectral densities for Nj = 1, 2 
and 4 flavours and have shown how they behave quantitatively differently when the sign 
problem is switched on. 
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We have also introduced a measure for the severity of the sign problem and have 
analysed the envelope, amplitude, number, and frequency of the oscillations of the spectral 
density once it becomes negative. 

It would be very interesting to compare the predictions we have made with ab initio 
lattice computations, in both settings of low and high density. There are some remarks 
to be made for this symmetry class, however. At zero chemical potential two-colour QCD 
has a substantial probability to have small or zero Dirac operator eigenvalues, which is in 
contrast to other QCD (-like) theories. Therefore it may be very difficult to reach small 
enough masses (that make the inversion of the Dirac determinant difficult) in order to see 
the effect of unquenching in the low-lying spectrum. 

The chemical potential introduces an extra repulsion of the spectral density from the 
origin, which might actually help in this situation. The sign problem on the other hand 
could be kept under control if the non-degeneracy of the quark masses is kept sufficiently 
small. 
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A The weakly non-Hermitian limit of the real eigenvalue density 
A.l Initial steps 

In this appendix, we consider the large- limit at weak non-Hermiticity of the real eigen- 
value density for the quenched case, which was not analytically determined in section 4.4 
of [34] . The real density of Wishart eigenvalues for finite matrix of "size" A^ (where A^ is 
even) and topology v is given by (see section 4.2 of [34]) 



Gn{x,x\h) 



(A.l) 



where 




(A.2) 




— oo 



dy sgn(x' - y) \x'y\-'^ e^-'^-'+v^ 2K,^{l^+\x'\) 2K^_{r^+\y\)lCN{x,y) 
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and the quenched kernel at finite N under the integral is given by 

N-2 , \ 2j / . , ,N, 



j=0 

For completeness we also give here the even skew-orthogonal polynomials since they appear 
in their limiting form in eq. (2.9), 

q2k(z) = (Vr/_)^'=(2A:)! ^2.(^) • (A.4) 
We now define the weakly non-Hermitian large- limit by 

A) - ^lim ^ ^ ^; . (A.5) 

The problem is that the limit in eq. (A.5) and the integral in eq. (A. 2) do not commute. 
Using the fact that 

/oo / r^' \ 

dy sgn(x' - y) f{y) =[ dy + 2 dy - dy\ f{y) (A.6) 
-00 yj-00 Jo Jo J 

we can correspondingly write Gw{x,x'; fi) as the sum of three parts, 

'; fi) + 2B{x, x; fi) - C{x, x; fi)} . (A.7) 
For A and B, it is possible to interchange the limit and the integral, and so we have simply 

, hy,(x') 
A{x,x;fi) = {-iY- — 7 ,„ / dyK.yj{x,y)hw{y) , (A.8) 



[sgn(x')]''/2 J- 

B{x,x'-fL) = r dyJC^{x,y)hM, (A.9) 



[sgn(x' 

where the scaled weight function and the weak kernel are given by eqs. (3.17) and (2.8), 
respectively, 

\x\ 



8/1 



1 



K-w{x,y) = ^^^^^2 j dss'^e "^f"^-'' {y/xJy+i{s^/x)Jy{s^) - ^Jyj^i{s^)Ju{s^/x)] 

The main problem is then reduced to determining C{x,x'; fi), where the non-commutativity 
issue remains. 

Let us next perform some fairly straightforward manipulation of C{x,x';fi). We can 
immediately take N ^ 00 for all the A^-dependent factors that lie outside both the integral 
and the sum. Therefore 

C{x, x'; A) = r^¥^ J™ ^^(^; A) ' (A-10) 

dZ'Kil'^ [sgn(a;')J'^/^ Af-s>oo 



-31 - 




(A.ll) 



(x o y) I . 

Since is finite, we can switch the integral and the sum. We will also split the two terms 
in the curly brackets, writing Dj\[{x; ft) = D^{x; fi) — D^{x; fi), where 

N-2 / , \ 2j 



with 
and 



with 



.0-,.,. . ( (il±^) (^) .... . 

(A.15) 



2N ' 



A. 2 Solution for limit of Dn{x;(i) 

Our plan is to try to evaluate the large- limits of eqs. (A.12) and (A. 14) in a manner 
similar to the way we calculated the weak kernel in section 4.4 of [34]. This involves taking 
the limit for large j and A, with j = tN and t £ [0, 1], of each of the factors inside the 
sums in eqs. (A.12) and (A. 14), and replacing the sums with integrals. So, in the present 
case, our primary task will be to determine appropriate large-A limits of P{j, N, jX) and 

It turns out that for D^{x; fi), this method works in a relatively straightforward man- 
ner. However, for Dj^{x;fi), we also find an extra contribution that does not involve an 
integral. This arises from a different large- A scaling of Q{j, N, fi). 

A. 2.1 Determination of limits of P{j, AT, p,) and D^{x; p,) 

Let us start with eq. (A. 13). First, we use [51, eq. (8.432.3)] to replace the modified Bessel 
function of the second kind with an integral, giving 

/ 1 -I- \ ''^^ fOO fOO 



r(^) V ^^^^ 



/* OO /* oo 
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with 

After writing the Laguerre function exphcitly as a polynomial, we interchange the order of 
integrations and perform the integral over y. The resulting sum can be recombined into a 
binomial as follows, 



T{^) jl ^_AiV Ji (u-l)(-+3)/2 / ^(^y+-+i ■ 

(A.18) 

We now wish to set j = tN, and take the large- limit whilst keeping t fixed. In fact, j 
should remain an integer, so we will set j = [tN], where the square brackets denote the 
integer part. 

In order to take the limit, it is necessary to split P{j, N, jl) into two parts, -Pi(j, A, jl), 
where the integral runs over 1 < ti < 1 + jc^/N, and P2{j, N, ji), where the integral is over 
the range 1 + ji'^/N < u < oo. For Pi{j,N,ji), we write u + 1 = 2 + 0{N^'^), and then 
make a change of variables to 

^ = 1 - (A.19) 

N{u — 1) 

with the result that (dropping the ©(A^^) which will vanish when the limit is taken) 

i+AV^ (_l)i2(^"i)/2Ar('^+i)/2 /-oo sJ(l + s)(^-l)/2 / ^2n 
du ■ ■ ■ = — -; / ds — — — — 1 



1 /i'^+i Jo (2 + 3)^+"+^ V 2As 

(A.20) 

It is important to note that the overall sign depends on whether j is an even or an odd 
integer. This implies that no single limit of Pi{j,N,jL) will exist. To evaluate the integral 
on the right-hand side of eq. (A.20), we change variables to n = j/s, leading to 

(A.21) 

In fact, we are only interested in the large- A limit of this integral, with j = [tN] and t 
fixed. As A — 7- c«, the integrand in the curly brackets converges to an integrable function 
in a manner that allows us to interchange the limit and the integral. Therefore 

POO />oo p ( f + l \ 

hm / ds--- = e'^'* / duu^''-'^/^ = e'^'* • (A.22) 

JV-^oo Jo Jo 2(^+i)/2 ^ ' 

j=m -"^ -"^ 

Combining everything, we have 

lim ^i(i^Ar,/l) ^ 2^+2 V^;lt('-^)/2e^'*. (A.23) 
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Note that we have been somewhat imprecise with our notation; there are actually two 
limits, depending on whether j steps through the even or the odd integers. 
We now turn to consider the limit of P2{j, N, fi), for which we can write 

^2' 



du ■ ■ ■ 



due 



u 



1 



N 



(A.24) 



where 0(2;) is the Heaviside step function. We can then use the Monotone Convergence 
Theorem to show that 



lim 

j=[tN] 



^ (U- l)('^+3)/2 



du 



exp 



r,2+ 1 



2fiH 
u-l 



u+2 



(A.25) 



where we used a simple change of variables s = (n + — 1) in the last step, and En{x) 
is the exponential integral defined in eq. (3.18). Our formula for the limit of P{j,N,fi) is 
therefore 



lim 



2^+2 y^e'^'* 



+ {-iy , (A.26) 



where it should be understood that there are different limits depending on whether j takes 
even or odd values. 

Let us now consider the limit of D^{x;/j,). Because of the two different limits of 
P{j, N, ft), technically we should now split the sum in eq. (A. 12) into odd and even parts, 
and treat each part separately. However, because of the way that everything combines, it 
is completely equivalent if we merely use the "average" limit over even and odd j, 

I Pij,N,fl)\ 2-+^ 



lim 

' j=m 



r(^ 



2 



(A.27) 



where it should be understood that this is the average limit over a small interval around t 
which includes (an infinite number of) both even and odd terms in j. 
We can now proceed as in [34], 



N-2 ^ 1 

lim > — = / 



7V->oo N 
j=0 

2N \ 



dt. 



lim 

j=[tN\ V ^ ^ 2N . 



exp[— 2/i^t] , 



J 



hm , . , 

j=[tN] ' 



(A.28) 

(A.29) 
(A.30) 



lim 



L 



u+l 



X 



1 



j=[tN] 



mil 



2N ' 



7V^™L N^+i V4iv 
2Vt^ 



( 

Ju+iiVxt) 



(A.31) 
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Hence 



lim^I)+(x;A) = W) £ dte~f^'H^^+'y^E._^{f,H)J,+,iV^t). (A.32) 



A. 2. 2 Determination of limits of Q{j, N, fl) and fi) 

For Q(j, A^, /i), we find that we expect the hmit 

lim (A.33) 

j=[tN] 

to exist. The key point is that there are two powers of N difference compared with the 
P{j, N, ft) case, which we can see as follows. Using the recurrence relation [51, eq. (8.971.4)] 

zL-^+Hz) = U + l){mz) - LJ+i(z)} + ^iz) (A.34) 

we can show that 



Q{j,N,fi)_^r fi 



N \ 2N 



2 



) [{3 + l){p{i.Nr^)-P{j + l.N,fL)] + uP{j,Nr^)\. (A.35) 



Since j + 1 oc at large A^, there are two powers of N difference between the leading 
orders of P and Q, and not one, as might naively have been expected. 

Leading order contribution: We use eqs. (A. 26) and (A.35) to show that 

lim ^%^ = (-1)^2-^+5 V^/it(-+i)/2e'^'S (A.36) 

j=[tN] 

where the sign depends on whether j is even or odd. Note that the other terms from 
eqs. (A. 26) and (A.35) are of smaller order, and so will vanish in the large-A limit. 

Now let us consider the limit of D^{x; fi). The presence of the N'^^'^ scaling for Q is 
the crucial difference compared with the P case, which scaled only with N'^. Although the 
N'^^'^ here can be matched with the N^^'^ prefactor in eq. (A. 14), this means that there is 
then no residual A^^^ to make dt when we take the large-A limit (see eq. (A. 28)), and so 
we will not get an integral. 

Because adjacent even and odd terms have opposite signs, however, we can use the 
fact that, for a continuous function f{t) on the interval < t < 1, 



i=o 

to take the large- A limit of the sum in eq. (A. 14). This results in a contribution of 



Regime 1 



lim D^{x] fi) = ^ fi [e"^"* Vt Ju{^t)\ = ,^ ' Jui^/x) . (A.38) 



t=o 
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Next-to-leading order contribution: There will also be a contribution from the next- 
to-leading order term of Q. Here, we can average even and odd terms (as we did for 
iV, fl) in eq. (A. 26)), and this average gives zero for the leading order. So the effective 
limit (over a small neighbourhood of t) is here (from eq. (A. 35)) 



N^oo iV'+l / \ dt I \ N^oo N" 

■ j=[tN\ I \j=m 



where we used 



Hence 



(2^)^+^^*^+^'* [E^ifLH) - E._^{(,H)] , (A.39) 

^ En{x) = - En-i{x) . (A.40) 
ax 



lirn^ D-^{x-fi) = ^ £ dte-^'H^+^l^ [E^il^h) - E._^{ft)} MV^t) 



N^oo 
Regime 2 



(A.41) 



We should add that lower orders will not contribute to the integral when we take the 
large- limit of the sum in eq. (A. 14). 

A. 2. 3 Final formula 

Pulling everything together and simplifying, we have 

32^7r [sgn[x')Yi^ [ /i T (^) Jq 



.42) 



^ El-, (p^t) J,.+i (y/xt) - fi^Vt (e^i-. (ph) - El-, (p^t)) JjVxt) 1. (A 

22 \2 ^^JJ 

We then substitute this (after an easy change of variables), together with eqs. (A. 8) and 
(A. 9), into eq. (A. 7) to give eq. (3.16). The final result, although not obtained in a 
mathematically rigorous manner, is well supported by numerical checks. Furthermore, 
we have also verified analytically that, on taking the Hermitian limit /i — ?■ 0, our result for 
the density shows complete agreement with [41]. 

B Decoupling of heavy flavours 

Looking at the RMT partition function (2.1) - and in fact for any QCD-like theory with 
or without chemical potential - it is expected that the microscopic spectral density would 
reduce to that for a fewer number of flavours when one or more of the quark masses were 
sent to infinity, 

p(f^)(C) pL'^'-'ho pL'''-'\c) • • • , (B.i) 



^We recall that the masses are purely imaginary. 
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We note that an even number of masses must be sent to infinity in the strong limit. The 
decouphng can easily be seen by numerical methods. 

We do not give a full proof. Instead we wish to illustrate the decoupling for 

^(A.,^2,C)^^^^n^^(A.,^l,C)(^^_ (B.3) 

To show this, we use that the weak kernel /C^ obeys 

/C^(i;,m^) = -/C^„(m^,■l;) ~ /^„(rn^)g^(t>) for - > 1 , (B.4) 

where 

^"("^'^ ^ e"i'^'/.+i(|m|) . (B.5) 
Thus, in the limit —mf ^> 1, eq. (3.20) yields 



X 



-fw{rh()qnjim^) 

= pl''^=''''\^;m2), (B.6) 

as seen from eq. (3.25). These steps can be generalised to higher Nf in a straightforward 
manner. 

The decoupling occurs in the sectors of real/imaginary eigenvalues as well, as seen in 
figures 7 and 8. For illustration, let us prove 

Pw ' (g;"ii,m2) > pw ' "'{t.;m2). (B.7) 

With eq. (B.4) it can be proven that 

Gw{x,x') '^^~°°> fw{x)Qw{x') . (B.8) 
Using eqs. (B.4) and (B.8) in eq. (3.24), we find, in the limit —ml S> 1, that 

iNf=2,m), . 

Pw [C,, mi, 1712) 



w vsj "T f f ^ 2\ t - 2\ 



-qw{rhl) 

= Pw' (?;m2) , (B.9) 

as seen from eq. (3.26). Moreover using eqs. (B.4) and (B.8) once again we can easily prove 
p(f^=^'«^)(e;m2) ^ for -rnl » 1. 
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C Sign-quenched partition functions 

In two-colour QCD the fermion determinant has no complex phase but only a sign, so we 
shall introduce a "sign-quenched" partition function, in which det[T>{fi) + rh] is replaced 
by I det[V{n) + as a direct counterpart of the phase-quenched partition function in 
three-colour QCD. The purpose of this appendix is to derive analytical expressions for 
such sign-quenched partition functions, first for finite N, and later in the limits of both 
weak and strong non-Hermiticity. The results in this appendix will be used in the main 
text for the analysis of the sign problem. 

C.l Main result for finite N 

Let us begin with a formula for the probability measure, with Zj = xj + iyj, 

Nf N „ N 



[TTlf - Zj) 

f=l k=l £ j=l f=l 

[N/2] 

X E ;7i7]73^n'^'(^2.-i-.2.)'^(y2n+i)---<5M, (C.l) 

n=0 i=l 

which is a minor modification of [34, eq. (2.9)]. Whilst being equal to eq. (3.3) for even 
N, this expression is more convenient for our current purpose. The weight function w{z) 
is defined in [34, eq. (2.12)], and below we only use that w{z) = 'w{z*) for z £ C, and 
w{z) = h{z) for z G M, with h{z) defined in eq. (3.5). The prefactor reads 

= (/x-independent factor) x (2^)-a^(a^+^)j^-^(^+'^-i)/2 ^ (c.2) 

and we refer to [34, eq. (3.46)] for more details. 

We also define a "partially sign-quenched" partition function, which we denote by 

^r^'^'^^'V; {-},{-}), as being identical to 4"^^^^^^^ in eq. (C.l) except that the 
fermion determinant is replaced by 

N 



n 



N'j, Nf 



Y[ \rhh\''\rh\ - 2j| rn'jimj - zj) 
h=i f=i 



(C.3) 



After these prerequisites, we can now state our main result at finite A^, 



Civ-fc \x^\ Z'^ {f^;{m}) 

(C.4) 

Below we give a proof of this formula for k = 1. The generalization to > 1 is then 
straightforward. 
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C.2 Proof: the fc = 1 case 

By definition, we have 

Nf N ^ N 

R[''j(\z) = U^'fU d'z, w{z,)\A^{{z})\ n IlK - ^^■) 

(/^;{W)/=i k=il j=i/=i 

[A^/2] n N 

^ E nUN - 2n)! H'^'^^^.-i - ^2.) '^(j/2n+i) • ■■5{yN)^5\z - z,) . (C.5) 

n=0 '' i=l e=l 

For 2n + 1 < i < N, 6'^{z — zi) in the last sum will yield 6{y) (with z = x + iy), hence 

Nf N „ N Nf 

<5(y)45«(x) = U^'fU d'z, w{z,)\Ar,{{z})\ [] HK 

(^;{"i})/=i k=il i=i/=i 

[N/2] 1 " ^ 

^ E ^!(iv - 2n)\ n '^^(^2^-1 - 4) ^(j/2«+i) • • • ^iVN) Yl ^^i^-^e) 

n=0 i=l e=2n+l 

Nf N ^ N Nf 

[Jm^Jjy d^Zk w{zk)\AN{{z})\l{l[{mj - zj) 



(nif - Zj) 



2^^N^'''\^J'^Am}) f=i k=i 



i=i /=i 



J2 nUN-2nV. H '^^(^2i-i - ^aJ '^(y2n+i) • • • S{yN) S^iz - zn) 

n=0 '' i=l 



Nf N 



w{x) JJ("T'/ — x) m'f 



/=i 

N-l „ / Af-1 \ N-1 Nf 

X [] / d^zfe \an-i{{z})\ n - ^^i n n( 

fc=i / \ e=i J j=i f=i 



W/2] ^ 

^ E ^i(jv_ i_2nV n '^'(^2^-l - 4)'5(y2n+i) • • • SiyN^i)S{y) . 

n=0 i=l 

(C.6) 



Therefore 



R 



C.3 The large-AT Hmit of weak non-Hermiticity 

Here we shall compute Z^^ s=m,^) ^^-^ {m}) = (| det[P(/x) + m]|) in the limit of weak non- 
Hermiticity defined in eq. (3.11), 

I? = 2Ni? , 
rhf = l^/Nnif , 
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Putting Nf = and x = rr? G M in eq. (C.7), we find 

Af-s>oo wym'^) ' 

In the weak limit, it follows from eqs. (3.10) and (3.15) that 

lim R\''j(=''^\m^) ~ i?!f^=°''')(m2) = -G^m^^m^) , (C.9) 

lim ~ . (C.IO) 

Af^oo w(m^) hyj{rh?) 

Prom eq. (C.2), we have 

r (^,,\-(N-l)(N-l+v)-{N-l){N+u-2)/2 
hm ~ , ttt;^ = (2/i) ^ 7?, ((-^-11) 

~ V2iV/x e^'^' = fL ei"'/^ . (C.12) 
Collecting all results, we finally obtain eq. (4.23), 

^^-ll^ll-)(A;m)~-Ae'^V^'^ """"^ 



/i^(m2) 

C.4 The large- AT limit of strong non-Hermit icity 



Below we show how to derive an analytical expression in the limit of strong non-Hermiticity 
for "^"''^^(//; mi, 7712) = (| det[P(/i) + mi] det[P(^) + m2]|), where N goes to infinity 

with n and mf fixed (mj = mj). 

Putting k = 2, Nf = 0, xi = mf, and X2 = m2 (with mf, m| G M) in eq. (C.4) and 
using w{x) = h{x) for a; G M, we find 

(iV, = ||2||,.), , C^_2 1 |mim2r (jV,=0,R). 2 .o^ 

(/x,mi,m2)~ \ml-ml\h{ml)Kmir^^'' {rn„m,) . (C.13) 

Next, let us recall that according to [34, 52] the real two-point function reads 

^2^~°''*^(^i'^2) = h{xi)h{x2)fN{xi,X2) , (C.14) 

where 

fN{xi,X2) = /CAr(xi,X2) Sgn(xi-2;2) 

+ dx dx h{x)h{x') sgn(xi — x) sgn(x2 — x) — sgn(x2 — x) sgn(xi — x) 

X 1Cn{xi,x)1Cn{x2-,x') — -'ICn{xi,X2)1Cn{x,x') . (C.15) 

We can plug eq. (C.14) into eq. (C.13) to get eq. (4.25), 

{Nf=\\2\\,u), , , . /s(mf,m|) 

mf — mi] 
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where fsixi,X2) is defined as the large- Umit of (xiX2)^''^/Af (a^i, 3^2), with no additional 
rescaling of parameters involved, and is explicitly given in eq. (4.26) for nonpositive argu- 
ments. In eq. (4.25) we omitted a /i-dependent prefactor on the right-hand side, as our 
primary interest in the strong limit is in the mass dependence of the quantities considered. 
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